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Abstract 



The purpose of this article is to give a pot-pourri of results on the 
mechanical and geometrical properties of curves and explicit solutions 
to problems on trajectories of particles under suitable force laws. The 
factor which unifies these rather disparate results is the ubiquity of 
the expressions in the title. We show how this explains in a unified 
fashion a plethora of properties of a large class of special curves. We 
also introduce two related ideas, the so-called d-transformation key res 
of a function and a new duality between trajectories for central force 
laws and those for parallel laws. The former will be our principal 
tool for obtaining many results by a reduction to cases which can be 
solved by elementary methods; the latter will allow us to move back 
and forth between results for central force laws and corresponding ones 
for parallel laws. 

Our treatment begins with an overview of the theory of the motion 
of a planet under a centripetal law, in particular one which varies as a 
power of the distance from the centre. We derive directly a criterion 
for a given orbit under a centripetal force to correspond to a power 
law, This is essentially a modern version of Newton's method. 

It is implicit in Newton's treatment of the two-body problem that 
a single orbit suffices to establish the force law and we give a short 
and elementary proof that this is true in a strong form — namely that 
knowledge of the affine curvature on an infinitesimal segment suffices. 

We then use the above criterion to give a detailed discussion of the 
Kasner- Arnold duality between power laws. This involves the first two 
expressions of our title and the ^-transformation mentioned above. 

We proceed to discuss two related themes in the differential geome- 
try of special curves — the three differential expressions of the title and 
a duality between curves which have natural representations as spirals 
and those which can be conveniently described by parametrisations of 
the form (F(t),f(t)), where F is a primitive of /. We show that the 
former are ubiquitous in the computation of geometrical or mechani- 
cal properties of such curves and display two families, the MacLaurin 
spirals and catenaries, which are particularly rich in these respects. 

In the final section, we give explicit parametrisations of trajectories 
under a parallel force field, emphasising again the case of power laws, 
and show that the differential expressions of the title and their stabil- 
ity properties with respect to the ^-transformation provide a unified 
approach. We close with a brief discussion of the theme of recreating 
curves from their curvature functions. 
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1 Introduction 



We start the technical part of this paper in section 2 below with an overview 
of one of the most fascinating topics in the history of mathematical physics: 
the relationship between Kepler's and Newton's laws of motion. This was the 
central theme of Newton's Principia. We are particularly interested in the 
following two statements about an object (called the planet) which moves 
around a second one (the sun) subject to a centripetal force (equivalently, in 
such a manner that Kepler's second law holds — equal areas are swept out in 
equal times): 

I. If the orbits are conic sections (more precisely ellipses, parabolas or hyper- 
bolas with the sun at a focus), then the force is inversely proportional to r 2 
where r is the distance to the sun; 

II. The converse statement: if the planet moves under an inverse square law, 
then the orbits are as above. 

In the first edition of his monumental treatise [Ne] Newton derived the 
statements of I (treating the three cases separately) and stated that the 
converse , i.e., II, holds. In later editions he added a proof of this converse, 
albeit one which is the subject of some controversy to this day. 

Newton's solution to what is often referred to as the Kepler problem is 
generally regarded as one of the key moments in the history of science and 
represented the culmination of the work of Tycho Brahe, who collected the 
data on the planetary motions, and Kepler, who distilled his three laws from 
this data. 

In addition, Newton developed a criterion which allowed him to deduce 
from the geometrical form of an observed orbit whether the planet was mov- 
ing under a power law and to determine which power was involved. He used 
this to investigate further exotic orbits which satisfy power laws not of the 
Kepler variety. Two of the more remarkable facts that he deduced were that 
orbits consisting of circles which pass through the sun and conic sections with 
the sun at the centre also arise from power laws (these cases will be dealt 
with below). 

It is implicit in Newton's work that the geometrical form of a single orbit 
suffices to determine the force law. We shall demonstrate that if we observe 
an orbit of a planet moving under a central field, then we can deduce the 
force law from its affine curvature. 

Our statements above use (as do all standard treatments of the dynamics 
of planetary motion) several hidden assumptions, since from the observation 
of a single orbit we can, in the most general situation, clearly only deduce in- 
formation about the force law at the points through which the planet passes. 
Basically, we are assuming, in addition to Kepler's area law (equivalently, 
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that the motion is determined by a central force emanating from the sun): 

a) that the force is independent of time; 

b) that the force depends only on the distance from its source (the sun) , i.e., 
it is invariant under rotation around the sun. 

c) that the force depends on the radius in a (real) analytic fashion. Hence, if 
we observe a non-trivial orbit (that is, one which is not a circle with centre 
at the sun) this will determine the nature of the analytic function on a non- 
degenerate interval and so in its entirety. (In our statements below we shall 
always assume that we do not have such a circular orbit in order to avoid this 
case — any force law which satisfies a), b) and c) clearly has circular orbits as 
special solutions.) 

Note that if we assume c) then it is not necessary to observe a complete 
orbit to derive the force law. Any non-trivial segment will suffice. 
One of our central results can be expressed succintly as follows: 

A) The force is proportional to a power of the distance from the central 
point (the "sun") if and only if the affine curvature K a e is also propor- 
tional to a power of this distance. In addition, the appropriate powers 
are related in a simple fashion, namely if oc ^ <^> Kgg oc ^j-. 

B) More generally, we have a force law K oc 0(r) if and only if K a g oc 0(r) / r 
where is a suitably smooth function (we use the proportionality sign 
rather than equality to avoid inserting constants). 

Curiously, we have been unable to find any relation between force laws 
and affine curvature in the literature (with one exception — see below). It is 
known that Newton used what in modern terms would be called curvature 
arguments but he can hardly have used affine curvature since the concept was 
only invented in the last century (see [Sc2] for a history of affine differential 
geometry). There are various articles on Newton's use of curvature (see, 
for example, [Br], [Co2]) but they clearly refer to the classical (Euclidean) 
curvature which was available to him as the inverse of the radius of the 
osculating circle. 

There is, however, one result on celestial mechanics which does invoke 
(implicitly) Euclidean curvature, namely the theorem of Hamilton [Ha] which 
states that the motion is Keplerian , i.e., corresponds to an inverse square law 

K oc — if and only if the hodograph is a circle , i.e., has constant curvature. 

(Recall that the hodograph of the motion is the curve traced out by the 
velocity vector). 

We can incorporate this into our scheme as follows: 
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C) The motion corresponds to a power law if and only if the Euclidean 
curvature K hod of the hodograph is proportional to a power of the dis- 
tance from the corresponding planetary position to the origin. More 
precisely, we have K oc r a if and only if /thod oc r~ a ~ 2 . 

The relationship between the force law and the two curvatures can be 
expressed succinctly in the formulae: 

1 K 1 K 

«hod-«afT OC — , Khod-K OC — , OC T. 

r 6 r 2 K aS 

In the course of our investigations we stumbled in a natural manner on 
the three-parameter family of functions 

/(f) =p(cos(d(t-t )))a 

which have the property that the corresponding orbits rf(9) = 1 are induced 
by power laws. Thecrucial fact is that the members of this family have the 
special property that the result of an application of any of the differential 
expressions of our title leads to a function which is proportional to a power 
of /; this fact will be the key to our treatment. More precisely, if 

f(t)=p(cQBd(t-t ))v 

then 

f 2 +f 2 = P 2 r 2d+2 , f+r = -P 2 {d-i)r 2d+ \ ff"-f 2 = - P 2 dr 2d+2 . 

Thus the curves of the form r d cos(d6) = 1 all describe planetary orbits 
which are induced by power laws (we omit the p and to since these corre- 
spond to the simple geometrical operations of dilation and rotation respec- 
tively). These curves were introduced by the renowned mathematician Colin 
MacLaurin, who was aware of precisely this fact. Further literature searches 
revealed that they had been investigated thoroughly in the period of classi- 
cal differential geometry (geometry of special curves). Their equations are 
usually written in the form r n = sin(n#) and they are known under the name 
of the MacLaurin spirals (sometimes sinusoidal spirals). A central thesis of 
our paper is that these all arise from the simple case /(f) = pcos(t — t ) 
(a straight line) under what we call the ^-transformation which associates 
to a function / a new one fd where fd(t) = (f(dt))*. The latter's stability 
properties with respect to the differential expressions of our title explain in 
a unified manner many of these results. 

The two standard treatises on special curves ([Go], [Lo]) devote exten- 
sive chapters to the development of some of the remarkable geometrical and 
mechanical properties of the MacLaurin spirals. 



6 



A further theme of our treatment is a duality between the curves discussed 
so far and a second family of curves which have analogous properties with 
respect to a force parallel to the y-axis, in particular one which is proportional 
to a power of the distance to the x-axis. 
Three motivating examples are 

1. The cycloid (t — cost, 1 + sint) which is a brachistochrone and tau- 
tochrone for a constant force parallel to the y-axis, and a trajectory for 



2. The circle (— cost, sint), which is a geodetic for the Poincare half-plane. 

3. the catenary y = cosh a;, which is the form taken on by a hanging chain. 

We remark that each of these curves can be parametrised in the form 
(F(t), f(t)) where F is a primitive of / which is in turn a function of the 
above form , i.e., p(cos(d(t — to)) 1 ^ (this is obvious in the first two examples, 
the third one is more subtle). 

The connection to the first class of curves is that there is a deep analogy 
between the properties of curves with parametrisations (F(t), f(t)) (where F 
is a primitive of /) and those of the spirals rf(9) = 1. Both have remarkable 
properties when / has the special form used in the definition of the MacLaurin 
spirals. It will be one of the main tasks of this article to display the reasons 
for these two facts. 

We shall also be interested in other families of curves associated with a 
force law — brachistochrones, tautochrones or isochrones, catenaries, elastica, 
geodetics for suitable metric tensors, and for light rays in media where the 
index of refraction is proportional to a power of the distance to a central 
point. For the sake of conciseness we will simply refer to these curves as 
trajectories. 

As a sample of these properties we mention that the MacLaurin spiral 
r d = cos(d9): 



a 



— law. 

y 




2. is an orbit for K oc 



1 




3. is a brachistochrone for K cx 



1 




4. is a catenary for K oc 



1 



■d+2 ' 



7 



(see [Go]). 

3) and 4) are examples of solutions of problems of the type: minimise 
J r a ds for a suitable index a. This suggests that the functions of the form 

p(cos(d(t — to)) 3 will supply solutions to many concrete problems of the 
calculus of variations, a fact which we shall verify and explain from our 
unified point of view. 

Our last theme is a description of the trajectories of particles moving under 
parallel-force laws. Recall that the trajectories of a particle moving under a 
given force law are the solution curves of the differential equations 

^ = /(*(*),!/(*)) ^ = 9«t),y(t)), 

where (f(x,y),g(x,y)) is the force at the point (x,y). The most inter- 
esting cases are, of course, that of a central force, i.e., one of the form 
/(r)(cos 9, sin 9) (where r and 9 are polar coordinates) which we discuss in 
section 2, or a parallel force of the form (0, f(y)). Of particular interest are 
the cases where / is a power function, i.e., K oc r a in the central case or 
K oc y a in the parallel one. Note that the family of trajectories is unchanged 
if we multiply the force field by a positive constant. Hence in the above 
cases, only the sign of the constant of proportionality affects the family of 
trajectories, not its absolute value. 

Of course, the famous result of Galilei (the trajectories for a constant 
parallel force are parabolas) is, together with those of Newton, one of the 
key results in the history of physics. 

One of Newton's less well-known discoveries is the fact that a Dido circle 
(i.e., one which is perpendicular to the x-axis) is a trajectory for a y~ 3 law. 
This will follow from our results; but we can go much further than Newton 
here. Thus, we obtain all trajectories for such a law and can show that these 
curves are the only circles which are trajectories for a y a law. 

The study of trajectories generated by force laws was an area of very 
active research in the first half of the previous century and is associated 
with Edward Kasner and his students. These investigators were interested in 
general properties of families of curves which arise as such trajectories, rather 
than in the explicit form of the families for concrete laws. We think that it 
is of some interest to document the fact that we can write down explicitly 
all the trajectories, using simple elementary functions (and an integration), 
in the case of a parallel power law. 

It is interesting that for the better-known case of a central power law, 
there are just three indices for which all of the trajectories can be described 
explicitly using elementary functions: the Kepler case K cx r -2 , Hooke's law 
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K oc r (where the orbits are conic sections with centre at the origin) and 
the Cotes' spirals (K oc r~ 3 ) which we discuss in section 2. In general, the 
trajectories for a given force law form a three-parameter or oo 3 family. For 
the general central power law the MacLaurin spirals mentioned above provide 
an explicit oo 2 family of trajectories, but in the general case the remaining 
ones can as far as we know only be described indirectly as far as we know 
(using functions which can be determined implicitly after a quadrature). In 
the Kepler case for instance, MacLaurin's family only picks up the parabolic 
orbits. 

The special case of rectilinear motion turns out, perhaps surprisingly, to 
be more intricate, and it is interesting to note that Newton, in his Principia, 
devoted a whole section to this case, which he regarded as a limiting case 
of the planar one (for a central force). In this case, it is, of course, not the 
geometrical form of the motion which is of interest but its direct description, 
i.e. formulae for the position as a function of time. Here the results are less 
satisfactory in the sense that we have to use not just elementary functions 
but also the inverse of such a function. 

In the final section we give descriptions of curves which satisfy the condi- 
tion that their curvature is proportional to a power of the distance from the 
x-axis. 

Many of the results of this article are, of course, known; our contribution 
has been to provide a unified approach. However, we do give explicit formulae 
and introduce special curves. 

As mentioned in the text, we were led to consider the questions below dur- 
ing the course of a cooperation with T. Russell and the late P.A. Samuelson 
to whom we owe thanks for many fruitful discussions. We would also like to 
thank Iain Fraser who read through and commented on an earlier version of 
the text (and, in particular, for his suggestingthe word "quadrality" as the 
appropriate substitute for "duality" for foursomes). 

We now turn to the technical part of this article. We begin with a survy 
of the Kepler problem. 

2 The Kepler problem 

2.1 Kepler's second law 

In this section we prove the basic fact that Kepler's second law (i.e., that 
the area swept out by the planet in a given interval of time is constant) is 
equivalent to the fact that the force is centripetal. This well-known fact was 
proved by Newton in [Ne], but since it is central to our approach we give a 
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proof in the spirit of what follows. 

We assume that the orbit of the planet has polar form r = f(9) for a 
smooth non-negative function /. Usually / will be strictly positive, but in 
at least one case we will allow it to have a zero (i.e. for the orbit to pass 
through the origin). Note that we are not assuming that / is periodic, i.e. 
that the orbits are closed. 

In this situation, the motion is determined when we know 9 as a function 
of time. The equations of motion are then 

x(t) = /(0(i))cos0(f), y(t) = f(9(t))sm9(t). 
The area swept out in the interval from t to t is 

A(t) = 1/2 / f{u) du. 
Je(to) 

The component of the acceleration perpendicular to the unit vector 

(cos 9(t), sin 9(t)) 

is the scalar product of the vector (— sin 9(t), cos 9(t)) with the second deriva- 
tive of (f(9(t))(cos9(t),sm9(t)) and an elementary calculation shows that 
this is 

2f(9(t))9 >2 (t) + f(9(t))9"(t). 

On the other hand, we can use the fundamental theorem of calculus and the 
chain rule to see that 

A"(t) = f{9(t))f{9(t))9'(tf + l -f{9(t))9"(t). 

Hence A"(i) is a multiple of the component of the acceleration and so the 
vanishing of A"(t) (which is just the analytical expression of Kepler's second 
law) is equivalent to the vanishing of the force component perpendicular to 
the vector from the sun to the planet. 

2.2 The inverse square law 

In order to prepare the reader for what follows, we consider briefly the Kepler 
problem in its original form, i.e., for the inverse square law. We shall show 
shortly that if we write the equation of an orbit in the rather unusual form 
rf(9) = 1, then the force is proportional to a power of the distance if and 
only if / satisfies a differential equation of the form 

f{9) + f"(9) = cf a (9) 
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for some constant c and index a. Then K oc r~ 2 ~ a . (The left hand side 
of this equation is, of course, the first of our differential expressions). Thus 
an inverse square law corresponds to the case where a = and this gives a 
clue to why the Kepler universe is particularly stable: the above equation 
is then linear, in fact it is simply f(9) + f"{9) = c. Of course this can be 
solved immediately, and the reader will see that the result is precisely that of 
Newton. Since the argument works in both directions, we have thus shown 
the equivalence of I and II above, i.e., completed the Newtonian programme 
with respect to the inverse square law. Of course, we still have to derive the 
above equation and this we will now do. 

2.3 The Newton-Somerville equation 

Suppose that we observe one planetary orbit, which we now write as a polar 
equation rf(9) = 1. The reason for using this form rather than the more 
usual one r — f{9) employed above is that this leads to a significant increase 
in transparency and clarity in the computations. In the Keplerian case, 
where the orbit is a conic section with the sun at a focus, the equation is — 
for a suitable choice of coordinate system — r(l + ecos#) = 1 where e is the 
eccentricity. In particular, < e < 1 corresponds to an elliptical orbit, e = 1 
to a parabolic, e > 1 to a hyperbolic. The fact that f{9) = 1 + ecos# is 
a much simpler function than its reciprocal is a further clue as to why our 
choice of equation is more natural in this context. 

Consider now the configuration consisting of the dilations of the single 
orbit and the rays going through the origin, i.e., the level curves of the 
functions (in polar coordinates) u = rf(9),v = 9. 

This is a so-called ^-configuration (see [Col]) and it is easy to calculate 
that the mapping (x, y) h-» (U, V) where 

2 

U = T -f{9f and V = g{9) 

is area-preserving and has our two systems as level curves, where g{9) is a 

primitive of „...„ . 

It follows that if the area condition is satisfied, then (up to constants) the 
time is given by t = g{9) and we can use this to get 9 as a function of t by 
inverting g. In the interesting examples (e.g. in the Kepler case), one can 
compute g explicitly, but not its inverse. Fortunately, as we shall see, for our 
purposes we shall not require this explicit representation for 9 as a function 
of time. 
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We have gone into this argument in some detail since it displays the con- 
nection with the concept of an S'-configuration. Of course, Kepler's area 
condition leads directly to this expression for t. 

We can now compute a simple formula for — and so (by the inverse 

do 

d6 

function theorem) for — in terms of / and its derivatives. In fact, as the 
reader will easily check, 

The motion is now (in Cartesian coordinates): 

_ cosfl(t) _ sinfl(t) 

x[t) -f(0(t)y y[t) -f(e(t)y 

where 9{t) = g^{t). 

We now have the machinery we need to compute the derivatives of (x(t),y(t)). 
We differentiate the vector function (x(t),y(t)) and use the chain rule to get 
the velocity vector as a function of time (in terms of / and its derivatives). 
The result is 

v(t) = -(sm9(t)f(9(t)) + cos9(t)f'(9(t)),-cos9(t)f(9(t)) + sm9(t)f(9(t)))- 
Similarly, the acceleration vector is 

a(t) = -(cos9(t),sm9(t))f(9(t))(f(9(t)) + f"(9(t))). 

We can immediately read off from the expression for the acceleration (which 
is, of course, proportional to the force) that we have a power law if and only 
if / satisfies the above equation 

f(9) + f"(9) = cr(9). 

(If we are only interested in the Kepler case, we can stop here since this 
closes the gap in our considerations above). This equation is equivalent to 
one which can be found in [Ne] — see [Cha] for a detailed discussion. Because 
of the central role that it will play in our considerations and of the relevance 
of the constants, we denote it by ns(c, a) and call it the Newton- Somerville 
equation (the first statement in a modern form of this criterion which we 
have been able to trace is in [So]). A recent reference is [Po] — c.f. equation 
(6.4) there. 
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2.4 A criterion for a power law 



If we differentiate this equation we get 

f(e) + f"'(9) = caf a - 1 (9)f(9) 
and so we can eliminate c by division to get 

m + ng _ m 
m + no) a f(o) 

f'(0)(f(0) + f"(0)) ' 
We can reverse this reasoning to deduce that we have a power law if and 
only if the derivative of the expression 

fir + n 
w + /") 

vanishes, in which case F oc r 13 , where j3 = — 2 — a and a is the (constant) 

value of 7(7+70" 

The equation for the existence of a power law is thus 

f\f+n(ff+f' 2 +fn+fr-(ff+fn(f'(f+n-f(f+n) = o- 

Summarising, the force satisfies a power law if and only if / is a solution of 

/(/'" + /') 

this ODE and the power p is then —2 — jjjj — . 

Due to the possibility of the functions in the denominator having zeroes, 
one should perhaps regard this equation more as a heuristic principle. Thus 
for a given / one tests (e.g. with Mathematica) whether it is a solution of the 
above equation. If this is the case, one computes the appropriate quotient 
and verifies it for constancy. Possible zeroes of /' or / + /" can then be 
investigated with ad hoc methods. 



2.5 Some applications 

Using this machinery we can instantly check whether a given orbit corre- 
sponds to a power law. We used it to check all of the examples in Newton. 
We then introduced parameters into the equations of such orbits and exper- 
imented to find combinations which produce further examples. We mention 
two simple ones: 
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Example We used the above equations to determine when a circular orbit 
derives from a power law by testing the case of a circle with centre at (a, 0) 
and radius 1. We found that this satisfies a power law if and only if a = 
or a = ±1. Hence a circular orbit corresponds to a power law if and only if 
the sun lies on the circle or at the centre. Both of these cases were discussed 
in Newton's Principia, but we have found no indication that he knew that 
these were the only ones. This fact is of some historical interest since one of 
the hypotheses which Kepler considered and rejected was that the orbit of 
Mars was circular, but with the sun displaced from the centre. 

Curiously, when we tried to extend this to the case of ellipses, i.e., to find 
out if there are other possible positions of the sun (other than the known 
ones at the centre or at a focus) for a power law, the computations turned 
out to be too complicated to be completed by Mathematica. 

Of the examples which we computed, we mention one which had interest- 
ing consequences. We tested the orbit 

f(9) = (a + bcos(d9)) 1/d 

for a power law. The above expression was cobbled together out of such 
examples as f{9) = 1 + cos 9, f{9) = (1 + costf)" 1 , f(9) = cos 9, f{9) = 
(cos9)~ 1 . The reason for the unusual nature of the dependency on d will 
become clear later. Apart from known or trivial cases, this provided us with 
two families of suitable functions, namely 

f d (9) = (cos(^)) 1 ^ 

and 

g d {9) = {l + cos{d9)f/ d 

The family 

f d (9) = (cos(d9)) 1 / d 

is well known in the classical theory of curves and, as mentioned in the 
introduction, will play a crucial role in our treatment. Its members satisfy the 
power law F oc r ~ 3+2d . The second family g d (9) = (1 + cos(d9)) l l d coincides 
essentially with the first one because of the simple identity (1 + cos d9) = 
2cos 2 (f). 

2.6 Two simple cases 

We return to our central ODE 

/ + /" = cf a . 
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with constants a and c. In general, this equation is non-linear. 

However, there are two cases where it is linear, namely a = (the Kepler 
case) and a = 1 (the K oc r -3 case), but it is linear for different reasons,this 
is evident from the different role of the constant, and it is instructive to 
compare the solutions. 

a = 0: the equation is / + /" = c. This is linear but inhomogeneous, and c 
plays the role of the inhomogeneous term. This equation can be solved by 
elementary methods and the solutions are exactly of the form 

f(0) = c + acos(6) + 6sin(#). 

(This completes our analysis of the Kepler case since these correspond to the 
polar equations of conic sections with focus at the origin). More precisely, 
if c > these / describe ellipses, parabolas or hyperbolas with the origin 
as a focus, and the freedom in the choice of the constants a and b (together 
with dilations) means that a curve specified by a solution of the differential 
equation, i.e., a curve of the form rf{6) = p where / is a solution of the 
differential equation and p is a scaling factor, is a curve of the required type. 
For the remaining values of c we get uniform motion in a straight line (c = 0), 
or hyperbolas (the Coulomb case c < 0). 

There is, however, one more case where this dependence can be computed, 
namely the case a — 1. Then the equation is / + /" = cf or /" = (1 — c)f. 
Once again this is linear, but c plays an entirely different role, namely in 
the coefficient of /. The solutions can be computed very easily and are 
logarithmic spirals, hyperbolic spirals and epispirals, depending on whether 
c < 1, c = 1 or c > 1. 

In this case we see that the class of orbits defined by the force law depends 
in an essential manner on the constant of proportionality. This case was 
settled by Cotes and the spirals which arise are known as Cotes' spirals. For 
a modern treatment, see Whittaker [Wh]. 

As a final remark in this section we mention briefly the more general case 
of a force law which is a sum of two powers (also of considerable historical 
interest). These will appear again in the next section. 

Example : Consider the equation 

f + f" = af a + bf 
which corresponds to a force law of the form 



■2+q r 2+/3 ' 
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The case a = and /3 — 1 is interesting both for physical reasons and also 
because it is the only genuine case of a sum of two powers which can be 
computed with ease. 

In this case the differential equation is / + /" = a + bf i.e. /" = (b — l)f + a. 
We distinguish the cases: 

a) b < 1. The solution is 

« + Ae^ 1/2e + Be-^ 1/2d ; 
l — o 

b) 6 > 1. The solution is 

+ Acos[(6 - l) l ' 2 6] + B sin[-(6 - l) l / 2 6}; 

l — o 

c) 6 = 1. The solution is 

- + e + B . 

3 The Affine Curvature of an Orbit 
Determines the Force Law 

It follows from the formulae given above that the statement of the title of 
this section holds. Due to its intrinsic interest, we go into this in more detail. 

3.1 Further geometric quantities associated with or- 
bits 

There are two further quantities which are determined by the geometry of the 
orbit and which turn out to be relevant — these are the curvature functions 
k and K h of the orbit and of the so-called hodograph, , i.e., the curve traced 
out by the velocity vector. The latter was introduced by Hamilton, who 
showed that the presence of an inverse square law is equivalent to the fact 
that the hodograph describes a circle, i.e., a curve with constant curvature. 
These curvatures can be calculated from the above equations using the stan- 
dard formulae from differential geometry for the curvature of a parametrised 
curve. Since one has to differentiate the equations of motion three times 
to achieve this task, one obtains potentially highly complicated expressions. 
However, a miracle takes place and they simplify to the tractable and signif- 
icant formulae: 

_ fV + f") J_ 

K (/ 2 + f 2 ) 3/2 ' Kh /+/"' 
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(These formulae can be computed by hand — we also checked them using 

Mathematica). The equation f + f" = — can be regarded as a quantitative 

version of Hamilton's characterisation of the Kepler case, since the circular 
form of the hodograph (i.e. the constancy of n h ) is equivalent to the validity 
of the equation / + /" = c for / (we are tacitly assuming that / + /" is 



Thus the expression f(6) + f"{6) which occurs in the central ODE is the 
radius of curvature of the hodograph. Since this expression and the related 
one f 2 + f' 2 will be of crucial importance below, it is of interest that they can 
be expressed in terms of the curvatures n and n h , together with /, explicitly. 
In fact, 



Thus we see how two of the expressions from our title arise in a natural 
way. We now consider the affine curvature of the orbit. 

3.2 Affine curvature 

Within the context of affine geometry, there are four concepts of curva- 
ture (affine curvature, equi-affine curvature, central affine curvature and cen- 
tral equi-affine curvature) depending on which type of geometry is involved. 
These are characterised by the choice of Lie group to define the geometry 
in question (in the spirit of Klein's Erlangen programme) — the affine group, 
the equi-affine group, the central affine group or the central equi-affine group. 
For completeness we recall briefly the definitions of these groups. 

The affine group is the six-parameter group of all affine transformations 
of the plane. (It will be convenient to use the classical terminology and refer 
to this as an oo 6 -group). 

The central affine group is the oo 4 subgroup of those affine transformations 
which leave a given point S invariant (S for "sun") 

The equi-affine group is the oo 5 group consisting of those affine transfor- 
mations which are area-preserving. 

The definition of the central equi-affine group (an oo 3 -group) should now 
be self-explanatory. 

Each of these groups is associated with an appropriate notion of curvature. 
In view of Kepler's second law it is natural to use the last group in dealing 
with orbital mechanics. In order to avoid the unwieldy terminology "central 
equi-affine curvature" we will refer to this simply as the affine curvature and 
denote it by K a g. 



positive). 



f + r 



■12 
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We can give our main result a more intuitive content by recalling that the 
afhne curvature has the following direct geometric interpretation. We denote 
by S the origin, by P a typical point on the curve and by P' a neighbouring 
point. We let the tangent to the curve at P meet the ray SP' at Q. Then 
the affine curvature at P is twice the limit, as P' tends to P, of the quotient 
of the area of the triangle P'PQ by the cube of that of SPP' (see [Sc] for 
details). 



3.3 A computational proof 

We shall begin with a purely computational proof of the relationship between 
the force law and the affine curvature, since this requires no knowledge of 
affine geometry apart from the formula for K a g. The formulae which we 
develop will also allow us to compute some simple illustrative examples. We 
will bring a conceptual proof at the end of the article. 
The required formula is 

/t aff = x A x/(x A x) 3 

where x(t) is a parametrisation of the curve and we use the Newtonian dots 

to indicate differentiation. The wedge product x A y of vectors x = (xi, x 2 ), 

y = (2/1,2/2) is the determinant x±y2 —x^yx of the corresponding 2x2 matrix 

(see [6] for the above formula). 

tp 1 1 / < n \ < n \\ (cos(0), sin(0)) . 

If we plug the parametrisation (x(B),y(B)) = — — into the for- 

j{9) 

mula for Kgg, a simple computation leads to the expression 

^= (f(0) + f"(9))f(9). 

The curvature of the hodograph is — — — , as can easily be com- 

f{0) + f"{0) 

puted by using the following parametrisation, which we include since it is of 
independent interest (it allows us to compute the hodograph of any planetary 
orbit from its geometric form): 

sin - f(9) cos 6, f(9) cos(^) - f'(9) sin(0)). 

We remark that this is the image of (f'(0), f{9)) under the reflection matrix 

cos 9 sin 9 
sin 9 — cos 9 

We remark that the MacLaurin spirals provide a large class of curves 
whose affine curvature is proportional to a power of the distance from the 
origin and which are therefore orbits for power laws. 
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Remarkably, they also have the property that the Euclidean curvature is 
proportional to a power of the distance to the origin — a fact that follows 
from the above relationship and the formula for the curvature of the orbit. 

Our results establish a relationship between the orbits of power laws and 
classes of curves characterised by the fact that the affine curvature is propor- 
tional to a a power of the distance to the origin. There are three cases where 
one can give simple explicit solutions to the first problem and hence to the 
second one. These are the Kepler case of an inverse square law, the Hooke 

case (K oc r) and the case where K oc — . As noted above, the second case is 

a well-known result, characterising the curves with constant affine curvature 
as the conies with centres at the origin. We have been unable to trace the 
other two in the literature and so state them here for the record. The conies 
with foci at the origin are characterised by the condition K a e oc r~ 3 . The 
third family are the so-called Cotes' spirals, which are characterised by the 
fact that their affine curvature is proportional to r -4 . They have the repre- 
sentation rf(9) = 1 in polar coordinates, where / is a solution of an equation 
f{9) + f"{9) = cf(9) for some c. The form of the solution depends on whether 
c> 1, c = 1 or c < 1. Typical examples are the spirals r cos(a#) = 1, r9 — 1, 
and logarithmic spirals. 

3.4 A conceptual proof 

We round off this section with a conceptual proof of the result on the affine 
curvature. As is to be expected this is very simple and follows more or less 
directly from the definition of the curvature notion that we use (which makes 
it all the more surprising that this result has not been documented already). 
Recall that the classical ,i.e., Euclidean, curvature of a curve is defined as 
follows. One introduces a special parametrisation, by arc-length, and shows 
that the second derivative of the parametrised curve is then proportional to 
the normal to the curve at the given point. The curvature is then defined 
to be the corresponding proportionality factor. In the case of the affine 
curvature one uses, as special parametrisation, the area spanned by the ray 
from the origin to a given point P on the curve. In view of Kepler's second 
law, this parametrisation is proportional to time in the case of a planetary 
orbit. The second derivative of the parametrisation is then parallel to the ray 
OP. The affine curvature is defined to be the corresponding proportionality 
factor (once again, the precise argument and formulae can be found in [6]- 
see §13). Since the second derivative is proportional to the acceleration, this 
explains the formulae used in our result. 
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4 The Kasner- Arnold duality 



4.1 Sinusoidal spirals and the ^-transformation 

We now explain the curious property of the sinusoidal spirals which we noted 
above, namely that they all arise from power laws. We begin with some 
remarks on these curves. They are usually specified by their polar equations 

r n = cos n6 or r n = sin n6 

where n is a parameter which can range over the real numbers, the second 
version explaining the nomenclature. The first version will be more conve- 
nient for our purposes (of course, the two are related by a rotation). 

This family was introduced by MacLaurin; particular choices of n produce 
some familiar classical curves. For example, n = 2 is a hyperbola (with 

centre at the origin), n = — 1 is a line, n = — - a parabola (focus at the 
origin), n = — - is Tschirnhausen's cubic (also called Catalan's trisectrix or 

1 1 

L'Hospital's cubic), n = — Cayley's sextic, n — - is the cardioid and n = 2 
the leminiscate. 

In a certain sense the family is generated by the two special cases n — ±1, 
i.e., the curves r = cos6> and r = (a unit circle tangential to the y- 

COS0 

axis at the origin and the line x — 1 respectively). Since this fact is of 
some moment in what follows we explain the process. If we employ the 
substitutions R = r n , = n9, then the equation R = cos© reduces to that 
of the generic sinusoidal spiral. Now the point with polar coordinates (R, 0) 
is the image of (r, 9) under the mapping z i— >■ z n of the complex plane. In 
other words, the sinusoidal spirals are the pre-images of the circle R = cos 
or the line R cos = 1 (depending on the sign of n) under this mapping 
(alternatively the images under z i— > z 1 ^ 1 ). 

At this point we remark that the correct setting for this discussion is 
not the punctured plane (i.e., R 2 or C without the origin) but its universal 
covering surface, i.e., the Riemann surface of the logarithmic function. This 
avoids the usual difficulties with the non-uniqueness or discontinuity of the 
argument function. However, in order to keep our dicussion elementary, we 
shall ignore this subtlety. 

Since we have written our curves in the form rf(0) = 1 (rather than 
r = g{6) which is more natural in the context of curve theory), we shall 
replace n by — d in the above equations, which now take the forms 

r(cos^) 1/d = 1 
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as before. 

These considerations lead naturally to the concept of the d-transformation 
fd of a function /, where 

/„(*) = (/(dt))*. 

As we shall see shortly, this transformation has particular stability properties 
with regard to the differential expressions of our title and this will explain 
many phenomena on the mechanical and geoemetrical properties of special 
curves. 

We now return to our central equation 

f(e) + f(6) = C f(er 

, which corresponds to a force law with j3 = —2 — a. Owing to its significance 
and also the importance of the constants a and c, we denote this equation 
by ns(c, a). 

4.2 The general situation 

The starting point of the ensuing discussion was a fact which struck us on 
reading Newton's Principia, in the version of Chandrasekhar [Ch]. One of the 
fascinating results he obtains is that when the orbits are conic sections, but 
with the Sun at the centre, then this also satisfies a power law, albeit with 
/3 = 1. These orbits are the images of the Kepler ones under the mapping 
z i — y z 2 (cf. Arnold [Ar]). This so-called duality, which apparently was first 
studied systematically by Kasner although we have been unable to find a 
precise reference, has been investigated in detail by many authors, see, for 
example, the same reference. The two families above suggest that there is a 
certain stability for power laws under transformations of the form z h> z n and 
it is this phenomenon which we now investigate. Simple computations show 
that this is not the case in general, but that for certain special situations, 
including the ones mentioned above, unexpected coincidences occur in the 
formulas and this explains the above facts. 

Thus our central question is: suppose that / satisfies ns(c, a). Does g also 
satisfy such an equation (with different constants) where g(9) = (f(d9)Y^ d , 
i.e., g is the d transform of /? Our above remarks indicate that there is a 
result of this sort for the special cases: 

a) f(9) = cos 9 (d arbitrary); 

b) f(9) = 1 + ecos9 (d = 2). 

As mentioned above, the expressions f{9) + f"(9) and f(9) 2 + f'(9) 2 will 
play an important role in our computations. 
In the above special cases their values are are: 
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a) f(6) + f" (9) = 0, /(fl) 2 + /W = l; 

b) + = 1, /W 2 + f{ef = 2/(0) + (e 2 - 1). 

Now there is a relationship between these two expressions. Indeed 

(This will be discussed in more detail below.) 
It follows that if / is a solution to ns(c, a), then 

a + 1 

for some constant b (or /(fl) 2 + f'(9) 2 = 2cln/(0) + b if a = -1). We 
emphasise that b depends on the particular orbit that we are examining, 

within a universe with power law K oc 2 and "gravitational constant" c. 

Those orbits with 6 = will be of particular interest below. For the case of 

a — law, these are precisely the parabolic orbits. 

A simple computation shows that if g(9) = f(dQ) l l d (and so f(d9) = 
g(9) d ), then, whenever / is a solution of ns(c, a), we have 

g\9) + g'\9) = -^-g^ d ~ d (9) + bg(9f~ 2d 
a + l 

for the constant b above and so 

9(0) + 9"(0) = C(2 + fr d) <7(0) 1W + 6(1 - d)g(9y-* d 
a + l 

for a ^ — 1. (The case a — — 1 must be dealt with separately). 

From these equations we can garner a wealth of information. Firstly we 
see that the dual of a power law is always the sum of two powers (again for 
«/-!)• | 

A particularly interesting case is when a = 1 and n = —. Then we see 

that a - law is dual to one of the form 

r 

2l + 2l 

as discussed above. 

There are three special situations where the dual collapses to a power law. 
a) c = 0, i.e., the case of straight line motion. Then g is a solution of 
ns((l — n)b, 1 — 2n) for any n. This explains the case a) above where f(9) = 
cos 9 and the dual curves are the sinusoidal spirals. 



22 



b) 6 = 0. Then g is a solution of 




c(2 + an — n) 



a + 1 



, 1 + an — n 



for any n. 

, 2 

This means that z \- > z transforms a — — 3 power law into a 2d — 3 

a 

law. This can be used to investigate the duality phenomenon in more detail. 
Since a great deal has been published on this topic we will not dwell on it 



the duality between conic sections with sun at the centre and such sections 
with the Sun at a focus, i.e., with (3 = 1 and j3 = — 2 respectively. Another 
interesting case is where n — — 1. This maps the case (5 = — 5 onto itself and 
thus shows that this case is self-dual (cf., Arnold [Ar]). 
As a final remark we note that the computations for the case a — — 1 are 
slightly different owing to the presence of the logarithmic term and we leave 
them to the reader. As a consequence we see that a dual in this case is never 
a power law (except in the trivial case d — 1). 

This is significant as it allows us to eliminate the third possibilty for 
dualities. 

c) where the two powers 1 + ad — d and 1 — 2d in the above equation coincide. 
However, this is the case where a — — 1 and so no non-trivial duality arises 
in this manner. 



5 A Duality between Trajectories for Central 
and Parallel Force Laws 

As we have seen, the members of the three-parameter family 



are solutions of ns(c, a) for various values of the parameters c and a, and so 
the curves of the form r d cos(d8) = 1 all describe planetary orbits which are 
induced by power laws (we omit the p and to since these correspond to the 
simple geometrical operations of dilation and rotation respectively). 

The main theme of this section is a duality between the curves discussed 
so far and a second family of curves which have analogous properties with 
respect to a force parallel to the y-axis, in particular one which is proportional 
to a power of the distance to the x-axis. 
We recall the three motivating examples 




f(t) = p(cos(d(t - t )) 



1/d 
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1. The cycloid (t — cost, 1 + sint), which is a brachistochrone and tau- 
tochrone for a constant force parallel to the y-axis, or a trajectory for 

a — law. 

2. The circle (— cost, sint), which is a geodetic for the Poincare half-plane. 

3. The catenary y = coshrr, which is the form taken on by a hanging 
chain. 

(We remark that each of these curves is a solution of a classical problem of 
the Calculus of Variations, for which see below. 2. and 3. are also trajectories 

for — laws). The connection to the first class of curves is that there is a deep 

analogy between the properties of curves with parametrisations of the form 
(F(t), /(£)) (where F is a primitive of /) and those of the spirals rf(9) = 1. 
Further, both have remarkable properties when / has the special form used 
in the definition of the MacLaurin spirals. It will be one of the main tasks 
of this article to display the reasons for these two facts. 
We can summarise the above remarks as follows: 

A) The curves of the first family are "spirals" with equations of the form 
rf(6) = 1 for suitable functions / of one variable. 

B) The curves of the second family have parametrisations of the form 
(F(t), f(t)) for a suitable function / with F a primitive of /. 

C) Further, when the functions / which occur have the form f(9) = 
(cos d9)d for some parameter d (or can be written in this form by 
simple transformations), then the curves have remarkable properties. 
As we shall see shortly, this is because these functions are the solu- 
tions of three particular ordinary differential equations involving the 
differential expressions of our title, a fact of some consequence for the 
mechanical properties of the corresponding curves. 

5.1 Curves of the form (F(t), f(t)) 

We therefore consider in more detail parametrised curves of the form 

(x(t),y(t)) = (F(t),f(t)) 

where F is a primitive of /. 



24 



In a certain sense, every generic plane curve with parametrisation (x(s),y(s)) 

( f* \ f s x'(u) 
can be reparametrised in the form / f,f\. We simply set t — / du 

\J J J y( u ) 

dt x'(s) 
which means that — = . 

ds y{s) 

Then if f(t) = y(s) and F(t) = x(s) we have 

= »w = /w- 

In particular, if our curve is the graph (s, y(s)) of a function, then we have 



f 1 , 

t — ; era. 



(«) 

We refer to this representation (which is essentially unique) as the canon- 
ical parametrisation of the curve. 

Note that there is no difficulty if we confine attention to curves which lie in 
the open upper or lower half-plane and whose velocity in the direction of the 
x-axis never vanishes. If the curve touches or crosses the x-axis or violates 
the second condition, then care is required and such a parametrisation may 
not exist. (Think of lines parallel to the y-axis). 

We illustrate this with a very simple example: the parabola (s,s 2 ). This 

has canonical parametrisation (— -, — ). Further examples will be computed 
below. 

Note that the parabola, which is originally "in one piece" , now splits into 
two parts (corresponding to t > and t < respectively). This mirrors the 
fact that the parabola touches the x-axis at its vertex. 

5.2 The duality (quadrality) 

These considerations make the following concept of duality between curves 
of the above two families natural: 

A) From spirals to parametrised curves: The spiral rf(6) = 1 corresponds 
to the parametrised curve I / f(u)du,f(t) 



B) From parametrised curve to spiral: If we are given the curve (x(s),y(s)) 
we compute the canonical parametrisation (F(t), f(t)) (i.e., f(t) = 
fx'iu) 

y(s(t)) with t = — — — - du) and then associate to it the spiral rf(6) = 
J U\ u ) 
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Of course, the transition from the spiral to the parametrised curve is 
immediate, while the reverse transition requires the intermediary step of 
computing the canonical parametrisation. We shall find that this procedure 
sometimes leads to rather surprising results. 

We remark at this point that it is an abuse of terminology to talk about 

a duality between curves since one normally does not distinguish between 

congruent or even similar curves. However, for our duality the position of the 

curve with respect to the system of coordinate axes is of crucial importance. 

Thus the dual of a unit circle with centre at the origin is very different from 

that of one which passes through the origin. This mirrors the fact that the 

x-axis and the origin respectively have a privileged role with respect to these 

two classes of curves. Further, the dual of a curve depends on whether we 

regard it as a parametrised curve or as a "spiral". Another aspect of the 

duality which will be of some consequence later is the fact that although the 

transformation from a curve to its dual does not arise from a transformation 

of space, nevertheless the transformation does act pointwise on the curves 

themselves. Specifically, the point (F(t), fit)) on the parametrised curve is 

. , (cost, sint)) . 
associated with —— on the dual spiral. 

/(*) 

For many purposes it will be useful to extend this duality to a quadrality 
by adding the inverse of rf(6) = 1 in the unit circle (i.e., the curve r = /(#)) 

and its dual, i.e., (^J g(u),g(t)^j where g(t) = y^y- Thus we have the 



scheme 



r/(0) = l (F(t)J(t)) 

t t 
r = f(0) ^ (l-}—du, ' 



m 7w. 

The main purpose of our note is to explain the significance of this duality, 
in particular for families which are trajectories in the above generalised sense, 
one with respect to physical properties with regard to r a laws, the other with 
respect to y a laws. 



5.3 Examples of duality 

I. We begin with duality for spirals. Note that we are using the term spirals 
in the rather loose sense of any curve with an equation of the form rf(8) = 1. 
Since the computation of the dual of a spiral requires only an integration, it 
is an easy task to set up a program, say in Mathematica, which automatically 
computes the inverse and both duals and displays plots of all four curves. We 
will therefore only mention a few cases which are of particular interest. We 
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begin with the Archimedes spiral r = a0, i.e., where f(9) = — . Its inverse 

ao 

is r = — , that is, f{6) = aO, which is known as the hyperbolic spiral. 
au 

The dual curves have parametrisations I — hit, — ] and I — -, at I respec- 



dualises to Neil's parabola ( "^">^ 2 



s a at J \at 2 

tively, i.e. they are the graph of y — -e~ ax and the parabola y 2 = 2ax. 

a 

II. As noted above, the canonical representation for the parabola (s, s 2 ) is 
( — \i J2) an d so the dual spiral is r = 9 2 , with inverse spiral r6 2 = 1 which 

III. The logarithmic spiral has the form rae w = 1, i.e., has f(9) = ae be . Its 

T 

inverse is -e~ w = 1 and so is again a logarithmic spiral. 

a 

Their duals are the curves with parametrisations (^ e& *' aebt ^j an d 

-—e~ bt , -e~ bt J , i.e., rays of the lines y = bx and y = —bx respectively 
ab a J 

IV. The generalised parabola y = ax a . In this case, the reparametrisation is 

(((l-a)t)i^,((l-a)t)T^). 

When a < 1 the parameter t ranges over the positive half-line, and when 
a > 1 over the negative half-line. 

This is a special case of the so-called higher-order spirals, which are written 

Q 

classically as r k = —. (For our purposes, k can be any real number, not 

2,71 

necessarily an integer). Its inverse is also a higher order spiral, and the dual 
curves are higher order parabolas, i.e., of the form y oc x a for a suitable a. 

Of particular interest is the case k = 2, which is known as the Fermat or 
parabolic spiral. 

Further interesting spirals are the Galilei spirals r = a — b6 2 , the logarith- 
mic spirals r = ae w and the rose curves r = a sin W. The special rose curve 
with b = 1 is one of the Cotes' spirals and is known as the epispiral. 

It is natural to generalise the rose curves to the family r d = a d sin(bd6) 
which generalises the family of MacLaurin spirals and has some of their 
propeties, albeit in a weaker form (for example, they are orbits for force laws 
which are the sum of two powers ). 
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5.4 The sinusoidal spirals and catenaries 

The specific properties of the sinusoidal spirals 

r d cos(d6) = 1 

suggest that their duals, i.e., the curves with parametrisation 

(F d (t)J d (t)) 

where fd(t) = (cosrft)s and Fd is a primitive of fd, will also possess interesting 
properties. It turns out that this is, in fact, the case, as we shall see shortly. 
We propose to call these curves MacLaurin or sinusoidal catenaries since 
they arise from the classical catenary using a transformation which will be 
discussed below. 

The key to this phenomenon lies in the fact that the functions f(t) = 
a(cos d(t — to))* have special properties — namely they are solutions of differ- 
ential equations of the form 

f{tf + f'{tf = af a 

f(t)+f"(t) = bf 

f(t)f"(t)-f'(t) 2 = cP 

as we shall see below. 

The expressions f 2 + f' 2 , f + /" and //" — f' 2 occur frequently in compu- 
tations involving the geometric and mechanical properties of curves (we have 
already met the first two). Hence those cases where they are proportional to 
a power of / (and so to a power of the distance to the origin or to the x-axis 
respectively) can be expected to have special properties. 

It is not a coincidence that the above family of functions satisfies each of 
these equations since the latter are closely related, as we shall see shortly. 

We note that the usual rule of thumb prepares us to anticipate oo 3 or oo 4 
resp. oo 4 solutions for these equations (one or two constants of integration 
and the two parameters in the equation). Thus we would expect our list to 
include all of the solutions in the first case, but not in the other two. In fact, 
it is also possible to give all of the solutions in the third case explicitly using 
elementary functions, and we shall do this in the next section. 
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6 On the expressions f 2 + f' 2 , f + f", ff" - f 2 , 
how to solve the equations f 2 + f' 2 = af a , 
f + f" = bf, ff - f 2 = cp and why one 
might want to 

We are interested in the three differential equations 

(AO P + r = af a ; 
(BO / + /" = bp; 

(co /r - r = cp. 

It will be convenient to consider the following more general situation: 

(A) p + r 2 = m 
(b) / + r = hp 

(C) ff - p 2 = p(f) 

where 0, and p are given smooth functions and the unknown / is a func- 
tion of one variable (so that our original equations correspond to the situation 
where the functions 0, ip and p are powers of / (with a constant of propor- 
tionality). 

The reason for this interest and the relevance of these equations for geo- 
metric and mechanical properties of our special types of curves will be dis- 
cussed below. 

The above stability property is a consequence of a simple computation 
which shows that if g is the ^-transform of /, i.e., g(t) = f(dt) 1 / d ) then 



fi~\p + n 

fi-\{p + n + d(fp - f 2 )) 
f*- 2 d(fp-n. 

Hence if / is a solution of (A), (B) or (C), then g is a solution of the cor- 
responding equations with the expressions g 2 ~ 2d <${g d ) , or g l ^ 2d {<p{g d )g d + 
dp(g d )) and dg 2 ~ 2d p(g d ) respectively as the right-hand sides. 



g 2 + g' 2 = 
g + g" = 
gg"-g' 2 = 
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6.1 The relationship between the equations 

These equations and their solutions are closely related, as we shall now see. 
1) If / is a solution of (A), i.e., if f 2 + f' 2 = 4>(f), then, differentiating and 
simplifying, we get 

2//' + 2/7" = 0'(/)/' 

and so / + /" = \4>'{f), i.e., / is a solution of (B) with ijj(f) = \4>'{f)- Also 
f' 2 (f)(f) 

1 + — = , and so, again by differentiating and simplifying, we get 
/ / 

(fr-f*) = \mf)-<Kf), 

i.e., / is a solution of (C) with p(f) = \f<t>'{f) - <f>(f)). 

Now assume this / is a solution of (B), i.e., f + f" = ip(f)- Then by the 

above f 2 + f' 2 = 20(/), where is a primitive of ip. Hence 

ff-f a = \mf)-4>u\ 

an equation of the form (C). Finally, if / is a solution of (C) then we can 
solve the (linear) differential equation 

\mr> - w) = tin 

for 0, given p, and so obtain an equation of type (A) for /. 

6.2 Some solutions 

In fact, using standard elementary techniques for solving O.D.E.s by quadra- 
ture, we can express the solutions implicitly as follows: 

(A) f 2 + f 2 = <f)(f) f(t) = F-\t + c), where F is a primitive of (0(w) 
(Ai) f 2 + f' 2 = af a special case of (A) 

(A 2 ) f 2 + f' 2 = a f a + d special case of (A) 

(B) / + f"=il>(f) special case of (A) with = 2^ 
(Bi) / + /" = bf special case of (B) 

(C) //" - f 2 = p(f) f(t) = exp(G- 1 (t + c)), i.e., In /(f) = G~ l (t + c) 

where G is a primitive of (2 J p(e u )e~ 2u )~^ 
(Ci) ff"-f' 2 = cp special case of (C). 
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We repeat for the record that if fit) = p(cosd(t — t ))d then 



j2 _|_ jl2 _ p 2 j-2d+2 

f+r = - P \d-\)r 2d+1 
ff-r = - P 2 df- 2d+2 . 

Thus this / solves (Ai) with a = p 2 , a = —2d + 2 
(Bi) with b = -p 2 , /3 = -2d + 1 
and (Ci) with c = -p 2 , 7 = -2d + 2. 

It is interesting that the hyperbolic trigonometric functions have similar but 
not identical properties, namely if g(t) = p(coshd(t — t ))d, then 



9 2 -g' 2 = 

9-9" = 

II 12 

99-9 = 



p 2 g- 2d+2 

-p 2 (d-l)g 

p 2 dg~ 2d+2 . 



-2<2+l 



We will mention a possible application of these facts below. 

The noteworthy form of the above solutions shows that they all arise from 
the simple function cost by applying the ^-transformation f(t) -w f(dt)d 
which, we recall, corresponds to the transformation z 1— > z d of the complex 
plane. This fact and, in particular, the fact that the latter is conformal, is 
relevant for some of the remarkable properties of the families (see, e.g. [Sc2] 
where this is used to explain many properties of the MacLaurin spirals). 

The reason for the above stability property lies in a simple computation 
which shows that if 
git) = f(dt) l / d , then 



9 2 + g' 2 = f*-\f 2 + f 2 ) 

9+9" = f l *- 2 ((f 2 +n+d(ff-n) 

99" -g' 2 = fi- 2 d{ff"-f 2 ). 

Hence if / is a solution of (A), (B) or (C), then g is a solution of g 2 + g' 2 = 
g 2 - 2d (f>(g d ), or g + g" = g l - 2d {ct>{g d ) + dp(g d )), or gg" - g' 2 = dg 2 ~ 2d p{g d ). 

In order to motivate these considerations we shall now show that the above 
three differential expressions are ubiquitous in calculations which are related 
to mechanical or geometric properties of curves. 
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6.3 Where the expressions arise 

As mentioned above we have two central motivations for these considerations. 
Firstly, the expressions of our title occur in many formulae which arise in an- 
alytically describing geometrical or mechanical properties of curves described 
as above. Secondly, and as a consequence, for curves which correspond to 
functions with the property that they are solutions of the above equations, 
then these quantities will take on a particularly simple form and so the curves 
will have remarkable geometric and mechanical properties. We believe that 
this is the natural explanation of the special properties of the MacLaurin 
spirals and catenaries. 

We bring a sample of the kind of occurence that we have in mind: 
For the spiral rf(9) = 1, then 

f 3 (f + f") 

1. the curvature is 

(f 2 + f'T* 

2. the afhne curvature is f 3 (f + /") (see below); 

3. if a planet moves around the curve in such a manner that it obeys 
Kepler's second law, then the absolute value of the acceleration vector 

is f 2 (f + n 

4. the curvature of the hodograph is (/ + 

5. the infinitesimal length is ds 2 = J^iP + f' 2 )d0 2 . 

For the curve with parametrisation (F(t), f(t)), 

ff" - f' 2 



1. the curvature is 



(/2 + / , 2) 3/2> 



2. if it moves under a force parallel to the y-axis which depends only on 
the distance to the x-axis, then the absolute value of its acceleration is 

ff" - f' 2 

p > 

3. The infinitesimal length is given by the formula ds 2 = (f 2 + f' 2 )dt 2 . 

A further situation where the expression f 2 + f' 2 occurs will now be discussed 
briefly. 
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6.4 The calculus of variations 



One of the most important aspects of the sinusoidal spirals and catenar- 
ies is that they provide a plethora of solutions to natural problems of the 
variational calculus. Elementary treatises on this subject typically use the 
following model applications, which are mathematically interesting, are of 
great historical interest and can be computed explicitly: the catenary (the 
curve which minimises potential energy and whose surface of revolution has 
a minimal-area property), the cycloid (brachistochrone), Dido circles (curves 
of a given length with endpoints on the x-axis which enclose maximal area). 
We shall now show that these are special cases of a general phenomenon 
and that this explains many of the remarkable properties of the MacLaurin 
spirals and catenaries. 

The common structure of these problems is that they maximise or min- 
imise functionals of the form 



J y 5 ds or j \ 



y ds or J r 7 ds. 

We consider the second case: If we use our basic parametrisation 

(cos(t), sin(t)) 

W) 

then it reduces to the variation problem of minimising or maximising a func- 
tional of the form 

J Fp(x,f(x),f'(x))dx 

where Fp(x, u, v) = ■u /3 (-u 2 +t> 2 )^ for a suitable (3 (f3 = 5 resp. (3 = —7—2). We 
now use the simple fact that for each j3, f is a solution of the corresponding 

Euler-Lagrange equation if and only if f 2 + f' 2 = af a for some a where 

0. 

f3 + 2 = — . This is seen by using the standard fact that in the case, as here, 

where the corresponding integrand is independent of the first variable, the 
Euler equation has the form: vf 3 — / is constant. This expression is 



I « 
CL2U 2 



when f 2 + f' 2 = af a , and so is constant when the above relationship holds. 

This explains many of the remarkable mechanical properties of the MacLau- 
rin spirals. 
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For example, in [Br] it is shown that the oo 2 family of parabolae of the 
form 

a(y — a) = (x — d) 2 

are the paths of light rays in a medium with refraction index oc y^. This is 
a special case of the above result since the canonical parametrisation of the 
above curve is 

( n t , ( t 
I 2a tan - + d, a I cos - 

and for f(t) = (cos |)^ 2 , we have / 2 + f' 2 = f 3 . 

We now show that the MacLaurin catenaries are also solutions of suitable 
variational problems. Interestingly, although the concrete variational prob- 
lems which are usually featured as text book examples are usually formulated 
for the case of a parallel force, many more concrete examples are known for 
the case of a central one. However, using the methods expounded here, we 
can extend the list of such examples for parallel forces considerably. 

We consider the problem of maximising 

J y 5 ds 

i.e. the variational problem with kernel 

J(y,y>) = y r (l + y> 2 )-^. 
The corresponding Euler equation reduces to 

y r {l+y' 2 )-h= c 
where c is an arbitrary constant. 

If we now use the parametrisation (F(t), /(£)), this reduces to f r+1 (f 2 + 

Hence if f 2 + f' 2 is proportional to f a , then the above parametrisation 
provides a solution whenever r + 1 = 2. Hence we obtain a solution with 
f(t) = p(cosrft)^ when 2 — 2d = a. 

6.5 Two more subtle examples 

As we have seen above, Dido circles and cycloids are examples of MacLaurin's 
catenaries. But why are catenaries, i.e., the graphs of the functions a cosh (-) 
also members of this family, as the name implies? This is because they have 
canonical parametrisations (F_i(i), f-i(t)) as the reader can verify. 
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A further fascinating class of curves are formed by the two parameter 
family of parabolae of the special form y = a + ^ x ~^ , which are the paths 
taken by light rays in suitable media. Once again, these are MacLaurin 
catenaries, explicitly (F-2(t), f-2(t)). This explains the following remarkable 
fact from [2]: The above family of parabolae are the paths followed by light 
rays in a medium where the index of refraction n(x, y) is proportional to y~? . 
Our method allows us to solve the corresponding problem for n(x, y) oc y a 
for any a. 

6.6 Miscellanea 

In this section we bring some miscellaneous remarks on our main topic. 

6.6.1 The explicit form of the parametrisation (F d (t) , f d (t)) 

The reader will have noticed that we never required the explicit form of 
the primitive of fd and so the explicit parametrisation of the MacLaurin 
catenaries. For the record we include it here: 
x(t) = - I ^(cos(d0 1+ ^cosec(^)F 1 2 (|(l + 1(3 + ±), (cos(rft) 2 ) 

y(t) = (cos^dt)) 1 * 

(we used Mathematica to obtain this expression — F% is the hypergeometric 
function.) 

6.6.2 The connection with contact geometry 

The representation of a curve in the form (F(t),f(t)), which might seem 
rather artifical and ad hoc, is natural within the framework of contact geom- 
etry as was pointed out to me by Valentin Lychagin. Namely, if we consider 
the simplest contact manifold R 3 = {(x, y,t) : x, y,t G R} with the natu- 
ral structure given by the differential form uj = dx — y dt, then the Legen- 
drian submanifolds are the curves of the form ((F(t), f(t),t) i.e. our special 
parametrisation displays the curve as the projection of such a manifold onto 
the (x, y) coordinates. 

6.6.3 Curves with prescribed curvature 

Many problems in applied mathematics can be subsumed under the general 
one: given a function f(x, y) on the plane or a subset thereof, we are required 
to find a parametrised curve such that at each point (x, y) on the curve the 
curvature k there is given by k = f(x,y). In other words we are looking for 
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solutions of the equation 

x(t)y(t) - x(t)y(t) 



(x(t) 2 +m) 



f(x{t),v(t))- 



It follows from above that the sinusoidal spirals are solutions of this problem 
for the case k oc r a for corresponding a. This fact, which was known to 
the classical differential geometers (cf. [Go], [Lo], has been rediscovered (for 
special cases) several times (cf. e.g. [Si]). Thus MacLaurin had shown that 
the radius of curvature of the spiral is given by the formula 

R = 



(n + l)r n - r 

Special cases of these reults can be found in the modern literature. All 
these results follow easily from the theory developed here. Analogously, the 

Maclaurin catenaries are solutions of related problems for — laws. 

All these results follow easily from the theory developed here. Analogously, 

the Maclaurin catenaries are solutions of related problems for — laws, as we 

shall show below. 



6.6.4 The rectification problem 

tp / / \ / \ (cost sin t) . . . , . 

if {x{tj,y{tj = ^— j is the canonical parametnsation of the curve 

rf(9) = 1 then ds 2 = j A {f 2 + f' 2 ) dt 2 . If (x(t),y(t)) = (F(t)J(t)) then 

ds 2 = (f 2 + f' 2 )dt 2 . For ^_Ld u , J-^ we have ds 2 = ±(f + f 12 ) dt 2 . 

Hence the expression for the infinitesimal length of the spiral rf(6) = 1 and 

for its "diagonal dual" (G(t), g(t)), where g = — and G is a primitive of g, 

coincide. This implies that the length between two points on the spiral and 
the corresponding ones on the diagonal dual are the same. This is the back- 
ground for many rectification results in the classical theory. For example, if 

r = aO, f(9) = — and so the dual curve is ( — lnt, — ) , i.e., the logarithmic 

' J w aQ \a ' a,tj' ' 

— -~ — ).u. 

the parabola y 2 = lax. In the words of [Lo] this means that the rectification 
problem for the Archimedean spiral is equivalent to the analogous problem 
for the parabola. The rectification problem played an important role in the 
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history of geometry, for which see the same reference. Since this question 
seems rather innocuous from a modern point of view, it is perhaps of interest 
to recall the historical perspective as described there: 

Thus we see that Archimedes had solved two of the three fun- 
damental questions for his spirals, namely the determination of 
its tangent and of the area swept out by a radial vector moving 
along the curve, but did not touch on the rectification problem. 
This important task could not be solved by his contemporaries 
or immediate successors, and would only be successfully attacked 
nearly two thousand years later, independently by Cavalieri, St. 
Vincentius, Roberval, Pascal and Fermat" (our translation). 

The method used here immediately givey a plethora of such results. We 
mention only one such, due to Fermat: the rectification problems for the 

spiral r k = — — and the generalised parabola z k+1 — (k + l)p k x, when a k = 

2,71 

2knp k ([10], p. 435), are equivalent. 

6.6.5 The Weingarten mapping of surfaces of revolution 

We consider the surface of revolution 

4>(u, v) = (f(v) cos-u, F(v), f{y) sinw) 

obtained by rotating the curve (0, F(v), f(u)) around the x-axis. Then if we 
compute the coefficients of the two fundamental forms, we find that 

E = f\ F = 0, G = f + f' 2 

and 

L = f/(f + /' 2 ) 1/2 , M = 0, N = (f 2 - ff")/(f + J' 2 ) 172 - 
Once again, we see the prevalence of the three expressions 

f + f\ f + f" and/f-/' 2 . 

As a consequence, if we rotate the MacLaurin catenary (Fd(t), fdif)) around 
the x-axis, then all of the quantities involved, i.e., E, G, L and iV (and, as 
a consequence, the principal curvature and Gaufiian curvature), are propor- 
tional to powers of the distance from the x-axis. (F and M vanish, of course). 

Without going into details we note that the Christoffel symbols for these 
surfaces (and hence the geodetic equations) also take on a very simple form. 
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If instead of using a classical Euclidean rotation we employ a hyperbolic 
one, i.e. the linear mapping with matrix 

cosh 9 sinh 9 
sinh 9 cosh 9 ' 

we obtain corresponding results for the Minkowski pseudo-metric (i.e., the 
one corresponding to the quadratic form x 2 + y 2 — z 2 ), which suggests that 
such surfaces could be of interest in the theory of relativity (this follows from 
the relationships satisfied by the hyperbolic trigonometric functions displayed 
above). 



6.6.6 Remark 

The characterisation of the functions p(cos(d(t — to)) 1 ^ as solutions of an 
equation of the form f 2 + f' 2 = af a is implicit in MacLaurin's treatment. 
He defines the spirals as curves which satisfy the condition 

d9 

tan V = r— , 
dr 

where V is the angle between the vector from the origin to the point (r, 9) 
on the curve and the tangent there. This condition is equivalent to the 
differential equation 

d9 r"- 1 

dr v / a 2n _ r 2n 

(cf. Gomes Teixeira [Go], p. 259, but note that there is a misprint in the 
equation there). If we put f(9) = then the above equation reduces to 

j 2 _|_ /'2 — a 2n j2n-2 



7 The trajectories for parallel power laws 

We now turn to a special case of the following problem: given a force field 
(in our examples on a subset of the plane), determine the family of all tra- 
jectories of particles moving under this field. An analogous problem is: given 
a function on a subset of the plane, determine the two-parameter family of 
curves which are such that at each point on the curve, the curvature corre- 
sponds to the value of this function. (Such problems arise, e.g., in theory of 
elasticae — cf. [Si]). We will give explicit formulae for the genuinely planar 
trajectories of particles which move under a force law of the form K oc y a . 
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In a second section we investigate free fall along a line parallel to the y-axis 
for such laws. In the final section we show that analogous methods can be 
used to solve the curvature problem, again for the case where the curvature 
is proportional to a power of the distance from the We are using the 

term "explicit" in the following precise sense. In the first case, we show that 
the trajectories can be parametrised by functions which are called elementary 
in normal usage (we include the hypergeometric functions in this category). 
In the two second cases, we have to allow the use of the inverses of such 
functions. 

More precisely, we shall show that the trajectories for a y Q -law are the 
curves with parametrisations (F d (t) , fd(t)) , where fd is a function of one of 
four particular forms for suitable parameters a, b and c, F d is a primitive of 
fd and d is a + 1, where a is the index of the power law. 

The four forms of fd are 



1 



adH 2 



for fixed d, this generates a two-parameter or oo 2 family of curves (the pa- 
rameters are a and an integration constant). The parameter t ranges over 
the punctured line, i.e. R\ {0}. (f d is, of course, proportional to t~d ) but 
we leave it in the above form for reasons which will soon become obvious). 



Ut) 



a{b - c ) 2 e a ^ dt 



^abdt _|_ fDacdty 2 

a three-parameter family. Here t ranges over the whole real line. 



Ut) = 



a(b- c fe a{b+c ^ dt \ d 

(pabt pacdt\ 2 



a three-parameter family. Once again, t ranges over the punctured line. 



f d (t) = (ac 2 sec 2 (acdtf) d , 

a three-parameter family (a, c and an integration constant), t ranges over 
the set of values for which acdt lies in the open interval ] — | , | [. 

The reader will have noticed that there is a problem with a - law, i.e. 

y 

where d = and we discuss this case later. 
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These are parametrisations (in the differential geometric sense) of the 
curves traced out by a particle moving under the corresponding force law — 
they do not describe its actual motion, i.e., t is not time. The parametrised 
motion is (u, f o F~" 1 (-u)), and cannot, in general, be written out explicitly 
using elementary functions but also requires the inverse of such a function. 

We think that it is of some interest to document the fact that we can 
write down explicitly all the trajectories, using simple elementary functions 
(and an integration), in the case of a parallel power law. We remark that 
these integrations can be carried out by Mathematica and the results can 
also be expressed in terms of elementary functions (if one admits this status 
to the hypergeometric functions). We have added the explicit forms at the 
end of the section. In contrast to the case of parallel laws, there are, for 
central power laws, just three cases in which all of the trajectories can be 
described explicitly using elementary functions: the Kepler case K oc r~ 2 
mentioned above, Hooke's law K oc r (where the orbits are conic sections 
with centre at the origin) and the so-called Cotes' spirals (K oc r~ 3 ). In 
general, the trajectories for a given force law form a three-parameter or oo 3 
family. For the general central power law MacLaurin produced an explicit oo 2 
family of trajectories (later baptised as the MacLaurin spirals — cf. [Col]), 
but the remaining ones can only be described indirectly in the general case 
as far as we know (using functions which can be determined implicitly after 
a quadrature). Thus in the Kepler case, MacLaurin's family only picks up 
the parabolic orbits. 

The special case of rectilinear motion turns out, perhaps surprisingly, to 
be more intricate, and it is interesting to note that Newton, in his Principia, 
devoted a whole section to this case, which he regarded as a limiting case 
of the planar one (for a central force). In this case, it is, of course, not the 
geometrical form of the motion which is of interest, but its direct description, 
i.e., formulae for the position as a function of time. Here the results are less 
satisfactory in the sense that we have to use not just elementary functions 
but also the inverse of such a function. 

For a brief introduction to the curvature problem and its relationship to 
the Euler elasticae, see Singer [Si] . A similar proviso applies to our solutions 
here. 

7.1 Planar motion under a parallel law 

We will now briefly discuss the method which leads to the above formulae. 
Once one has them, then it is a routine application of the chain rule to show 
that the curves they describe actually have the claimed property, i.e., that 
they are trajectories for suitable laws. The only interesting fact is to show 
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how these formulae arise, and this we shall now do. Since we are considering 
y a laws, then it is natural to confine attention to curves in the upper half- 
plane. In the Galilean situation, we must consider four cases, depending on 
the position of the parabola with respect to the x-axis, exemplified by the 
graphs of the following functions: 

y = x 2 

(a full parabola touching the x-axis); 

y = x(l — x) 

with x between and 1 (an arc of a parabola with vertex pointing upward); 

y = x(x- 1) 

with two unbounded branches corresponding to x < and x > 1; 

y = x 2 + 1 

(a parabola which is disjoint from the x-axis). 

Of course, only 6) is physical, i.e., corresponds to an attractive force to- 
wards the The other three describe trajectories under a repulsive 
force (anti-gravity). 

We remark firstly that in this section we are excluding the cases of free fall 
parallel to the y-axis, which implies that we can parametrise the trajectories 
in the form (F(t),f(t)), where F is a primitive of /, as shown above. Our 
proof is a combination of the following simple observations: 

A) If we parametrise a curve in the special form (F(t), fit)) as above (i.e. 
where F is a primitive of /) and a particle moves along this curve under 
a parallel force law, then the acceleration at the point (F(t), f(t)) is 

L f{t)f"{t)-f{tf \ 
V Pit) )■ 

This is a simple consequence of the chain rule. 
Hence 
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B) We have a y a law if and only if / satisfies an equation of the form 
f(t)f"(t) - f\tf = bf{t) fj for constants b and p. The relationship 
between this /3 and the exponent a of the power law is very simple — 
(3 = a + 3. Of course, the force is attractive or repulsive (with reference 
to the x-axis) according to whether b is negative or positive. 

C) If / satisfies an equation as in B), then so does its rf-transform fd (with 
distinct parameters b and P). This follows from the equation: 

hti-f d 2 = dr 2+ Hff-n 

which we gave above. 

The crucial point is now the following rather unorthodox parametrisations 
of the parabola: 

D) The general parabola (more precisely, the part or parts above the x- 
axis) has canonical parametrisation (F (£),/(£)) where / has one of 
following four forms: 




(F(t), fit)) is then a parametrisation of the parabola y = a(s — b) 2 where we 
choose F so that F(0) = b. (Again, t ranges over the punctured reals and 
the vertex is "lost"); 
or 



= g{b-cfe aib+ ^ 

^ ^ ' ^ e abt _|_ e act^ 2 ' 

(F(t), fit)) is then a parametrisation of the parabola y = a(x — b)(c — x), 
ormore precisely of the arc in the upper half plane, t ranges over the real 
line; 
or 



_ a{b - c ) V( 6+c )* 
Jv-J ~ 



(F(t), f{t)) then parametrises y = a(x — b)(x — c), or more precisely the two 
unbounded branches in the upper half plane. They correspond to t < and 
t > respectively; 
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or 



fit) = ac 2 sec 2 (at) 

where (F(t), f(t)) parametrises y = a((x — b) 2 + c 2 ). t ranges over the open 
interval 

(The four different forms correspond to the four possible positions of the 

parabola with regard to the x-axis, as described above). 

In these parametrisations, a is assumed to be positive and b < c. 

E) in each of the cases in D) the function / satisfies the condition 

fit) fit) - fit) 2 = 2afit)\ 

We have thus verified Galilei's result that the parabolas are orbits under a 
constant parallel force law, albeit in a rather roundabout manner. However, 
the payoff from this approach now follows. 

F) If we combine this fact and the formula in C) then we see that the func- 
tions quoted at the beginning of the article (which arise from the above 
parametrisations of the parabola via the process considered there) sat- 
isfy the condition that //" — f' 2 is proportional to a power of /, in fact 
to f 2+d , and so the curves mentioned in the first paragraph are trajec- 
tories for K oc y d_1 . Of course, as mentioned above, one can compute 
this directly; but the treatment described here shows how to find the 
appropriate form for the functions /. 

We now obtain the results above as follows. We suppose that we have a 
trajectory for a y a law parametrised in the form (G(t), g(t)). Then for any a 
we can find a suitable d so that g(t) = fd(t) = f(dt) 1 ^ d , where / corresponds 
to the Galilean case, i.e., a constant force. We then use the above formula 
for the possible forms for / and so find g as above. 

It thus only remains to consider the singular case, that of a - law. In 

y 

this case, the differential equation for / (i.e. //" — f' 2 = cf 2 ) can be solved 
directly, and one sees that the trajectories are (F(t),f(t)), where / is a 
function of the form expp(t) with p a quadratic function. 

For the sake of completeness, we add the explicit forms of the functions 
Fd, i.e., the primitives of the corresponding fd for cases 2), 3) and 4) at the 
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head of the article (the case 1) is trivial): 



F d (t) 



F d (t) 




) 



(ac 2 sec 2 (adt)) d sm(2adt) 



These expressions were computed with Mathematica (in the first equation 
we have assumed that a = 1, b = and c = 1 to make the formula more 
tractable). 

In the singular case (i.e. a - law), the formula for F(t) requires the use 

y 

of the error function from statistics. Since the explicit form depends on the 
nature of the quadratic function, i.e., on whether it is definite or indefinite, 
we omit the details. 

As final remarks, we note firstly that despite the complex form of these 
formulae, they do, of course, specialise to simple curves for certain choices of 
the parameters. For example a Dido circle (that is, a circle with centre on 

the x-axis) is a trajectory for a — law (interestingly, this case was treated 

y 3 

by Newton in his Principia). Secondly, using the theory developed here, it is 
easy to derive a criterion for a given curve with parametrisation (x(s),y(s)) 



A(s) = (3y'(s)y"(s) + y(s)y"\s)) x'{sf - (3x» (s)y> (s) 2 + y( S )x'"( S )y'(s) 

+3y(s)x"(s)y"(s)) x'(s) + y(s)y'(s)x"(s) 2 
B(s) = x>(s)y>(s)(x'(s)y"(s)-x"(s)y>(s)). 



The index of the power law is then this constant minus three. 

Using this criterion one can check that not only are Dido circles trajectories 
for the y~ 3 power law, but also that, as mentioned above, they are the only 
circles which are trajectories for any power law. Ellipses and hyperbolae (not 
necessarily rectangular) with the x-axis as axis are also trajectories for y~ 3 
laws. 



to be the trajectory of a power law. It is that 



A(s) 



be constant, where 



B(s) 
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7.2 Rectilinear motion 



We now consider the case of rectilinear motion. This was treated by Newton 
in Section VIII of his Principia. Interestingly, he did not deal with it directly 
but considered it as a limiting case of planar motion under a central force. 
We shall do the same here, with the difference that we shall employ a parallel 
force. 

Firstly we note that this is equivalent to solving a differential equation of 
the form 

V" = F{y) 

for the specal case where F is a power function. Now there is a standard 
method for solving equations of the above type by quadrature: we introduce 

the variable w where y' = w(y). Then y" = w——, so the equation reduces to 

d 

— w 2 = 2F(y) and so we have w 2 = 2Fi(u) du, i.e., w = a/2Fi(m) du where 
ay 

F 1 is a primitive of F. 

Thus we have reduced to the equation y' = ^j2Fi(u) du, which we can 
solve by quadrature. (This is essentially the modern version of Newton's 
treatment as exponded by Chandrasekhar [Ch]). Using the methods of the 
first section, we can give a more direct approach which does not require 
quadratures. 

If we suppose that / is a function which satisfies the condition 

//" " f 2 = af a 

and we set y(t) — f o F _1 (£), (with F again a primitive of /), then the 
computations above show that 

d 2 y , 

and the same argument as we used above shows that we obtain all of the 
solutions in this manner. Once again, we can choose / as one of the special 
cases above and thus obtain analytic expressions for rectilinear motion under 
a parallel power law. 

The case of a - law is again an exception. 

y 
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